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Abstract 

We consider two-dimensional large gauge theory with D adjoint scalars 
on a torus, which is obtained from a D+2 dimensional pure Yang-Mills theory 
on T^"*"^ with D small radii. The two dimensional model has various phases 
characterized by the holonomy of the gauge field around non-contractible 
cycles of the 2-torus. We determine the phase boundaries and derive the order 
of the phase transitions using a method, developed in an earlier work (hep- 
th/0910.4526), which is nonperturbative in the 'tHooft coupling and uses a 
1/D expansion. We embed our phase diagram in the more extensive phase 
structure of the D + 2 dimensional Yang-Mills theory and match with the 
picture of a cascade of phase transitions found earlier in lattice calculations 
(hep-lat/0710.0098). We also propose a dual gravity system based on a 
Scherk-Schwarz compactification of a D2 brane wrapped on a 3-torus and 
find a phase structure which is similar to the phase diagram found in the 
gauge theory calculation. 
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1 Introduction and Summary 

Gauge theories on spaces with compact directions have been studied for a 
long time. As a prototypical example, d + 1 dimensional Yang-Mills the- 
ory at a finite temperature T corresponds to a compactification of the (Eu- 
clidean) time direction on a circle of length /3 = 1/T. It is obviously im- 
portant to study such a compactification to understand the physics of con- 
finement /deconfinement transitions p. More generally, one can consider 
Yang-Mills theory on a compact space S. If the volume of S is finite, there 
is no sharp phase transition, but for an SU{N) gauge theory in the large 
limit there are sharply demarcated phases depending on the shape and size 
parameters of S. In case the compact space is a torus, the phase diagram 
as a function of various radii (and coupling) reveals a rich phase structure 
including a cascade of phase transitions in which the "Polyakov" loops 
along various non-contractible cycles become non-zero in succession as the 
radii are reduced [H[5]. Most of these studies are numerical (in the lattice or 
in the continuum) or, in some cases, based on holography (see Section |6] for 
references and more details). One of the motivations of the present paper is 
to investigate these questions analytically in a simple situation, as explained 
below, by using and extending the "large D" technique developed in [B]. 

To elaborate further, let us consider a Euclidean rf + D-dimensional gauge 
theory on a (i + dimensional torus with radii H 



S 



dt 



n 

M=l 



dx 



M 



-TiF. 



d+D 



(1) 



Here the length of the temporal circle is denoted as Lq = /3 and the rest are 
denoted as Lm, M = 1, ...,d + D — 1. The phases of ([T]) are characterized by 
Wilson lines around the d + D noncontractible cycles of the torus: 



W., = TiU, 



2 / Af,dx^ 





(2) 



^In this paper, we will not consider the contribution of the 6 term. 

^Our notation for spacetimc coordinates is: {a;° = t, x^^}, M = 1, ...,d + D — 1. We will 
further split the d + D — 1 coordinates into d 'large' dimensions {xo,x'^},i ~ 1, ...,d — 1 
and D 'small' dimensions x^ ,1 ~ 1, 2, ...D (the meaning of 'large' and 'small' is explained 
below) . 
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where no sum over is intended. These Wilson loops transform nontrivially 
under the centre symmetrjl^l. For sufficiently large radii L^, all PV^ vanish, 
signifying unbroken centre symmetry. In this phase, local gauge-invariant 
observables are independent of in the strict large limit [9l [TOl fTT] . 
Since Wq can be interpreted as exp[— 6*5], where Sq is the action for a static 
quark, the phase with (Wq) = exhibits confinement. As is well-known, as 
/3 is reduced (i.e. the temperature is increased), below a certain critical value 
Pc, (Wq) becomes non-zero, signaling a deconfinement transition together 
with a breaking of the centre symmetry Zp'^ — )■ Zp'^~^. In this phase, 
the observables can depend on /3 but are still independent of Li [10]. It has 
been argued from lattice studies (see [1] and Section [6] for a review) that 
as the other radii are successively reduced, one has a cascade of analogous 
symmetry breaking transitions Zp~^~^ — )■ Z'^^~'^ — t- ... — )■ 1. 

While it would be fascinating to study all the above phases analytically, 
in this paper we will be able to study the phases oi & D + 2 dimensional pure 
Yang Mills theory on T^"*"^ (i.e. ([1]) with c? = 2) in which a D-dimensional 
torus (with radii / (27r), J = 1, 2, D) is taken as small (ensuring broken 
Z^v symmetries in those directions), leaving the remaining d = 2 directions 
(including, time) of variable size. Such a theory is given by a Kaluza-Klein 
reduction □ of ([1]) on the small T^, and is described by the following action: 

-J2^\Y',Y']lY',Y']]. (3) 

I,J / 

Here comes from the gauge field components and the covariant 

derivative is defined as = — i[An, ]. A naive KK reduction leads to 

^ 'Centre symmetry' [71 [S] is generated by quasiperiodic 'gauge transformations' a{x^) = 
exp[2ni{nf,xt'/Lf,)A], where A = diag[l/A^, 1/iV, ...,l/N, {1~N)/N]. The quasiperiodicity 
is up to phases = exp[2'!Tin^/N], fi — 0,1, d+D — 1 which parametrize d+D copies of 
the centre of SU{N), The a{x^') are vaUd gauge transformations locahy, and leave 

local colour-singlets e.g. tT^F^^ invariant; in particular they commute with the hamiltonian. 
However, under the a-transformations — )■ h^jW^j,. A non-zero value of (M^^) implies 
spontaneous symmetry breaking of the centre symmetry in the ^-direction. 

^ Kaluza Klein reduction is tricky for gauge theories [3J [TT] , since in the confined phase 
the KK modes can have energies ~ 1/{NL), which become arbitrarily low at large N. 
The fractional modes, equivalent to the 'long string' modes of [12], can be understood as 
arising from mode shifts of charged fields in the presence of Wilson lines whose eigenvalues 
are uniformly distributed along a circle (see Section [2] for an explicit verification for this 
statement). In the deconfined phase, however, the KK modes have energies like in 

ordinary field theories, and KK reduction proceeds as usual. 
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massless y^'s and g = g'; however, a mass m for the adjoint scalars as well 
as radiative splitting between g and g' is induced from loops of KK modes 
(see Appendix |X] and Section O respectively, for more details). 

We should remark that ([3]) can either be regarded as a step towards 
understanding the full phase diagram of the d + D dimensional theory ([1]), or 
be understood as a two-dimensional gauge theory in its own right. The spirit 
of the latter approach is to provide an example of an analytically solvable 
low-dimensional gauge theory in the limit of a large number of adjoint scalars 
(the d = 1 theory was discussed in [6]). Equation ([3]), from this viewpoint, 
provides a bosonic counterpart of the Gross-Neveu model [13] where D plays 
the role of Nj, and the S'0(D)-invariant bilinear Yl'i=i^a^b Section [271] 
plays the role of SU (A^j)-invariant fermion bilinears such as Ylfji ipiipi of the 
Gross-Neveu model. 

Note that in this second point of view, where we regard the action ([3]) as 
an independent theory in its own right, the mass m can be taken to be ar- 
bitrary. In particular, if the mass is sufficiently large (m^ ^ g'^ND, g''^ND), 
the phase structure can be determined perturbatively [HE]- However, as we 
will see in appendix if we regard ([3]) as a KK reduction of ([T]) on T^, then 
the mass m of the adjoint scalars is much lower than the scales mentioned 
above and the theory is not amenable to such perturbative methods. One 
of the goals of the present work is to provide a nonperturbative ^ analysis 
of ([3]) valid for any value of mass (including m = 0), based on the 'large 
method developed in the previous work [B]. 

A few additional comments are in order: 

(i) KK gauge theories have important applications to phenomenology 
[El [m [TTl [18] . Theories such as ([3]) provide important toy models in this 
context. In particular, issues such as different running of gauge couplings 
in the compactified and decompactified theories can be examined in such 
models. We will encounter some of these issues in Section [2j 

(ii) Gauge theories on compact spaces can sometimes have gravity duals. 
The deconfinement transition in TV = 4 super Yang Mills theory on 5*^ x 5*^ 
[T^ [2U] is a weak coupling continuation of the gravitational Hawking-Page 
transition [2T1 [22]. Similar correspondences for two-dimensional supersym- 
metric gauge theories on tori were analyzed in [23t [3]FI . In this paper, we will 
obtain ([3]) with D = 8 from a Scherk-Schwarz compactification of a three 
dimensional super Yang Mills theory with sixteen supercharges, which cor- 
responds to the world volume theory of D2 branes. The latter theory has 

^in the 'tHooft coupling. 

^For an extensive list of correspondences between low dimensional gauge theories and 
gravitational systems, see [Ml [^ . 



4 



an AdS/CFT dual [21], which leads to the construction of a gravity dual for 
([3]) in a sense defined in Section As we will see, the gravity analysis 
will complement our knowledge of the phase structure from the gauge theory 
analysis. 

(iii) Large two dimensional gauge theories themselves are interest- 
ing objects in the context of string theory and QCD. For instance, confine- 
ment /deconfinement type transitions [2B] have been analytically found in 2D 
QCD with heavy adjoint scalars [H]. In addition, stringy excitations and 
glueball spectra have been obtained in 2D models in |29l [301 ED ESI [33] . 
Thus, these models are good laboratories for real QCD. Our study, in fact, 
has a direct relevance for |32]; we hope to return to the issue of glueball 
spectrum discussed in this reference. 

The principal result in this paper is the determination of some parts of 
the phase diagram of the two-dimensional theory ([3]) at weak coupling. The 
result is summarized in figure HI The second result is the gravity analysis 
which complements the phase diagram at strong coupling, which is presented 
in figure O 

The plan of the paper is as follows. 

In Sections [21 El and m we analyse the model ([3]) at weak coupling by using 
the 1/D expansion [M] developed in [6]. We find that the nature and the 
order of the confinement /deconfinement type transition depends on whether 
the size L of the spatial circle is large or small. For large L (which corresponds 
to the TiV = phase), we find (see Section [2]) a single first order transition, 
thus providing analytical evidence for earlier lattice studies (see Section 16.11 
for further details). On the other hand for small L (corresponding to the 
TtV 7^ phase), the analysis in [6] is valid and, as detailed in Section [3l the 
transition consists of two higher order phase transitions. The phase diagram 
is summarised in Section [H in figure [H and is in agreement with those from 
the lattice studies of [21 [5l [1]. We compare our results with these lattice 
studies and with [3] in Section [61 

To supplement our 1/D analysis of the gauge theory, we consider in Sec- 
tion [5l a dual gravity theory, obtained from D2 branes wrapped on a 3-torus 
with a Scherk Schwarz circle. Although the gravity results pertain to thermo- 
dynamics in the strongly coupled regime, we observe that the phase structure 
is qualitatively similar to that of the weak coupling gauge theory. This al- 
lows us to arrive at a conjectured phase diagram in figure [5l which suggests 
a particular way of connecting the phase boundaries of figure [H 

^One of the motivations for this work was to construct a gauge theory dual to a dynam- 
ical Gregory-Laflamme transition. We will discuss this issue in a forthcoming publication 
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In Section we comment on the dependence of the order of phase tran- 
sition on the topology of the compact space. In Appendix |X] we discuss the 
masses of the adjoint scalars which appear from integration of the KK modes 
at the one loop level. In Appendix |B] we fill in some details needed in Sec- 
tion [2] for integrating out the adjoint scalars. In Appendix [C] we discuss the 
influence of the mass term in on the phase diagram. In Appendix [D] we 
provide important details of our gravity analysis. 

2 Confinement / deconfinement type transition 
in large radius torus 

In this section, we will analyse the confinement /deconfinement type tran- 
sition in (j3]) for large L. First, (in Section 12. ip we will integrate out the 
adjoint scalars in a. 1/ D expansion, in a manner similar to [6], leading to 
the effective hamiltonian fl23p in the large L limit. Next (in Section fI72\i we 
use , at L — oo, the results of [HI ESI ESj, who studied this effective action 
in slightly different contexts, to determine the phase structure of our theory 
([3]) at large L. The justification for extrapolating their result to finite L, as 
detailed below, comes from the phenomenon of large volume independence 
[HI [TOl [11] which is valid as long as L is large enough to ensure Tr^ = 0. 

To keep the analysis simple, in this section we consider with m = 0, 
and defer the case of nonzero mass to appendix O As we will find there, the 
inclusion of the mass term does not change the qualitative structure of the 
different phases. 

2.1 Large D saddle point 

In this section, we will generalize the analysis of the + 1-dimensional gauge 
theory [6] to the d = 2 model (|3]) and show that if we consider the number 
D of adjoint scalars to be large, the theory can be considered to be in the 
vicinity of a large D saddle point and various quantities such as the free 
energy and the mass gap etc. can be computed around the saddle point in a 
1/D expansion. 

As mentioned above, in this section we will consider (|3]) with m = 0. 
Introducing an auxiliary field Bab, the path-integral of the gauge theory can 

^Although the saddle point arises in a manner similar to that in four-fermion models 
such as Gross-Neveu or Nambu Jona-Lasinio, the saddle point is complex. See [B] for 
details. 



6 



be rewritten as 



S{B, A, Y) = l^dt £dx ^^F,\' + i {D,Y^) ' - i\B,,YlYi + -1^B,,M~,]^,B,, 

^ (4) 

where A/" is a constant factor and we have used the following matrix, 



ab,cd 



■^|Tr[Aa, \c][K Ad] + (a O 6) + (c ^ rf) + (a ^ 6, c ^ (5) 



A" (a = 1, ■ ■ ■ , A^^ — 1) being the generators of SU (iV). Our approach, similar 
to the one-dimensional case, will be as follows. We will integrate out the y-'^'s 
to obtain an effective action for and Bab-, and find a saddle point solution 
for Bab (for given A^) in a large D limit. The effective action for will 
be essentially obtained by substituting the saddle point value of Bab in this 
effective action. 

As in [6], it is convenient to decompose Bab as the sum of a trace piece 
(independent of x,t) and an orthogonal part: 

Bab{t)=l^^5ab + g'hab{t,x), (6) 

where hab satisfies J dt J dx baa = 0. Such a decomposition, into a large diag- 
onal piece and a small off-diagonal fluctuation, will be a posteriori justified 
by finding a saddle point for Bat of the form (Bab) = i^o^^ab, where Aq is a 
real constant (depending only on the 'tHooft coupling). 
With this decomposition, the action reduces to 

S,nt = ~ J^^ dt dx (^-^babYX^y (7) 

where we have used M~j^^^5cd = ^ab/'^N [6]. Let us now take a large D limit, 
^,^'^0, N,D^ 00 s.t. ~X = g'^DN, \' = g'^DN fixed. (8) 
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To the leading order of this expansion, we can ignore the interaction term 
5*1,1*0 In that case, the integration of bat will contribute just a numerical 
factor and we will ignore it. Integrating over the y^'s in 5*1, we then get the 
following leading result for the partition function 




S,ff[A, A] = -^^^ + dt f dx^^F^,' + 5^e//, (9) 
oQ Jo Jo 

where SSeff is the 1-loop contribution from the F^-integration: 

SS,ff[A,A] = ^logdet i-Dl + A') . (10) 

It is difficult to evaluate this last quantity in general. However, using argu- 
ments similar to [3], we will find that: 

(a) If ^ A = g'^DN (this assumption will be justified in (fT9|l ). terms 
in 6Seff involving derivatives of A^, e.g. terms involving the gauge field 
strength and its covariant derivatives are suppressed. 

(b) If LA > 10 we can approximately treat Aq and Ai as commuting 
matrices. 

The argument for assertion (a) can be sketched briefiy as follows. Con- 
sider the simplest of the 1-loop diagrams (figure [1]) which contribute to ffTOj) 
and has only two external gauge field insertions: 



p 




Figure 1: A simple Feynman diagram contributing to ([TUj) . 



^In the one dimensional model {d ~ 1), the next order of the ^jD expansion has been 
evaluated in [S]. There, such l/L* corrections do not change the nature of the phase 
structure. We can expect that the same thing will happen in our two dimensional gauge 
theory also. Hence we do not evaluate the 1/-D corrections in this article. 

"'^''This assumption will be justified later in what we will define as the 'large L regime' 
L> Lc (see ([55)1 and (P^). since icA will turn out be large, using p^ . 
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For large A, the F^-propagator in the loop carrying momentum p can 
be expanded in powers of p^/A^. The first term in that expansion goes as 
ND{k^ki, — k'^g^y)/^'^ (in the Feynman gauge), where the factors of and 
D come from the F^-loop, and k is the external momentum going into the 
loop. This term amounts to a correction, to the Fqi/q"^ term, of the form 
(1 + 0(A/A^)), as claimed above. We will ensure below that A/A^ ^ 1. 

The argument for assertion (b) will follow a posteriori after we proceed 
with the assumption that are constant commuting matrices. Under this 
assumption, as detailed in appendix [Bl the one-loop term becomes 



6S, 



D(3L 
87r2 



Ar2 (_7rA2 + n {A^ + A^) log (A^ + A^) - ttA^ log A^) 



y |Tr([/V')|^^ J^^ 



(Av/(/3A;)2 + {Llf 
ill 



Here we have used a momentum cut-off A to regulate the determinant!"], 
and used the notation U = e^^^° and V = e^^^^. Ki is the modified 
Bessel function of the second kind. Using the asymptotic expansion Ki{z) = 
\/ t: I2ze~^ + ■ ■ ■ for large z (justified below) and omitting some irrelevant 
divergent terms, we get 



DN'^(3L 



eff 



Stt 
D 



(A^+A^)iog(i+A;)+A^iog(|;^ 



L./^e-^^|Trt/|V/3W^e-^^|TrV^|' 1 + 



(12) 



where the ■ ■ ■ terms represent higher order terms in e~^^, e~^^ . Since we are 
interested in large L, it is obvious why higher order terms in e~^^ should be 
ignored. We now have an a posteriori justification for having ignored terms 
involving commutators [f/, V]\ the smallest gauge invariant such term would 
have at least 2 f/s and 2 V^s and hence are expected to be of the same order 
as the U'^V'^ term in fill I) , and hence can be ignored [3]. Higher order terms 
in e~^^ are ignored because they will be small in the region of parameter 
space in which interesting phase structures will appear. We will ensure this 
at the end of Section 12.21 (see comment (a) below 



^^The cut off A should be smaller than Mkk, which is the inverse length scale of the 
_D-dimensional compactification torus. 
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In the L — )■ oo limit we can ignore the term in f[T^ involving V. Let us 
integrate A from under this assumption and derive the effective action 
for the [/-variable. The saddle point equation with respect to A^ reads as 




1 - 



2A(3 



(13) 



In the TtU = phase (which is realized for sufficiently large A 
determined implicitly from the equation 



A, 



27rA2 



log 1 



A2 



IS 



(14) 



This equation can be viewed as a renormalization condition which assigns a 
A-dependence to A' such that the physical mass scale Aq is held fixed. Here 
Aq plays a role analogous to Aqcd and its choice specifies the theory. We 
will, in fact, choose it in such a way as to ensure the conditional 



A ~ A' at A = M, 



KK- 



(15) 



As P decreases, TtU eventually becomes non-zero. However, near criti- 
cality (3 A 2> 1 (this is justified because of ( 1T9|) and (l32l) ). Therefore, we can 
solve (fT3|) for A in the form Aq + 0(exp[— /3Ao]). We will write the explicit 
solution only for A ^ A: 



A = An 



2ttAI/\' + 1 







/ /5 





(16) 



where 



An = \I^W 
Zn 



27rA2 



(17) 



This follows from the fact that the distinction between g and g' in ^ vanishes in the 
original pure YM theory ([T|). In order to fix the precise coefficient of the renormalization 
condition ([T5|). we need to evaluate the contribution of the mass ([5^ and the running 
of g and g' for scales larger than Mkk- The value of mass and the running of the 
couplings depend on the details of the higher dimensional theory. If we are interested in 
the two dimensional gauge theory ([3]) itself as in the comment (iii) in the introduction, 
the renormalization condition that replaces (|15l) is arbitrary. However, even in that case, 
the qualitative nature of the phase structures in this paper does not change as long as 
(A^ + TO^)/A » 1 is satisfied. 
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and W{z) is the Lambert's W function defined imphcitly by the equation z 
W{z)e^^'^\ Note that, by using an expansion W{z) = log z — log (log z) + - ■ 
for large z, we obtain 



An 




Imposing the renormalization condition ( IT5|) . we then obtain a relatioio 

> 1. (19) 



A^/A ~ ^ log 



A 



Thus we confirmed that the assertion (a) above is satisfied. Substituting 
f|T6|) in Seff in (|9]) and ignoring the term involving Tr^, we get the following 
effective action for A..: 



S{A) 

Dm 



C(A',Ac 



/3 

dt I dx 



2XN 



2 

01 




C(A',Ac 



/3LA 



57r 



ttAI 
A' 




where the terms ■ ■ ■ are higher order in the same sense as in ( ]T2|) . We have 
used j dxdt\TiU\'^ = L/3|Tr[/p which is correct up to derivatives of U which 
occur at 0(A/A^) and are small, as argued above. 

In the limit of L — )■ we can choose the gauge Ai = 0. Solving 
the Gauss's law condition in this gauge, we get Aq = {g^/dl)g where g = 
— ^[Y^,dtY^] is the charge density; in the large D saddle point, especially 
for large enough Aq, the condensate is static and temporal fiuctuations of g, 
and consequently, of Aq, are suppressed. As a result, we can write 



dx dt TrF^-^ 



f3 / dxTiid^AoY = f3-^ / dxTT\d^U\ 



The inequality follows from the condition A/M^^ <C 1, which is necessary for a 
KK reduction. This condition implies that the dimensionless 'tHooft coupling is small: 
Ad+2M^^' « 1. 

^"^ Although we are apparently taking L — > co here, as we discuss at the end of the 
next subsection, we can extend these results to finite L which is large enough to ensure 
vanishing ot TiV. Note that if TrV ^ 0, the 1/D expansion does not work at L — !• oo as 
argued in [37] , where such a situation is discussed in the presence of i?-symmetry chemical 
potential. 
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where 



U{x) = Pexp[i / dtAo{x,t)] = exp[il3Ao{x)]. 
Jo 



(22) 



By rescaling x ^ x' = \(3x we eventually get an effective action in terms of 
U{x) 



S/DN^ = C{X\ Ao) 



dx 



^Tr(|9.f/n-i,|Trf/r 



where 



27rA2/33 



-Ao/3 



(23) 



(24) 



Note that is a monotonically decreasing function of (3. 

Equations similar to ([2SD, (1241) have earlier been derived in [2] who con- 
sider a two-dimensional gauge theory with heavy adjoint scalars of mass m. 
In their equations m appears in place of Aq. We could have, in fact, derived 
the above effective action (123|) as follows: the large D saddle point generates 
a dynamical mass Aq for the adjoint scalars which turns the theory into a 
massive adjoint scalar QCD; once this is established, we can use the method 
of [H] with m — Aq, to arrive at fl23p . The agreement with the results of |14j 
provides an additional check on our derivation. Adjoint scalar QCD with 
large scalar mass has also been considered by [3] who have independently 
derived equations analogous to (123|) . flM|) : our method of derivation follows 
their derivation closely, except that our mass is dynamically generated, as 
mentioned above. The discussion in appendix O involving arbitrary m relates 
the two extreme cases of large mass and zero mass. 



2.2 The phase transition at Large L 

Phase transitions in the system fl23|) have been discussed in [TU [35], [36] . We 
will adopt their result to infer about phase transitions in our two-dimensional 
gauge theory ^ at large L (the range of L is defined in f[55i) ). For complete- 
ness, we will briefly review some of the results in these papers. 

If we regard the coordinate x in (j23l) as time, it becomes the quantum me- 
chanics of a single unitary matrix, a subject that has been extensively studied 
[M], [33- IS]- Phases of such a model can be described by the behaviour of 
the eigenvalue density 

1 ^ 

p{o.x) = -Y^5{e-e,{x)), (25) 

i=l 
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where exp {i9i{x)) are the eigenvalues of 

The hamihonian of this system (regarding x as time) can be written as 



H = J de ^-pv' + -p'j - e - ^, (26) 

where we have ignored the constant term C in (123|) . Here v = dgll, and 
n(^,a;) is the canonical conjugate of p{9,x). We have also used the notation 
I-'JyU"' , u_n = (wn)*- Note that Un{x) are moments of the eigenvalue 



u 



n 



density 

M^'^) = ^ E «n(x)e— ^ (27) 



27r 

n=— oo 



To study the various equilibrium phases of the system, we study static solu- 
tions of (126|) . which [HI |35l [36] are given by v{9) = and 

p{9) = ^ (^^/eT2^^^^) , (28) 

where pi = ui = ul is the first moment (real in this case), which must 
self-consistently satisfy (see ([2j 



/ d9p{e) COS 9 = Pi. (29) 
The constant E is fixed by solving the normalization condition 

/■2n 

/ d9p{9) = l. (30) 

Depending on the value of the constant ^ in (126|) . Eqn. (130|) may not deter- 
mine E uniquely. In general, we obtain three branches E{C,), depending on 
whether (I28p describes a uniform, non-uniform or gapped eigenvalue distri- 
bution (see figure [2]) • 
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— TV TT — TT TT — TT TT 



Figure 2: Configurations of eigenvalue density p{6) in tlie unitary matrix 
model. Tlie left plot is the uniform distribution (corresponding to pi = in 
fl28|) ). the middle one is the non-uniform distribution {\E/{2$,pi)\ > 1) and 
the right one is the gapped distribution (|_E/(2,^pi)| < 1). 

The value of E{C,) in these three cases are different. The thermodynamic 
stability for the various branches of the function E{C,) is analyzed [I ^ l35 | 136] 
by comparing the values of the Euclidean Hamiltonian (126|) which can also 
be regarded as the free energy. They can be summarised as follows: 

• Independent of the value of ^, the uniform solution always exists. We 
call this phase as I. 

• For ^ < ^0 = 0.227, only one solution (phase I) exists and is stable. 

• At ^ = .^0) there is nucleation of two gapped solutions. One is unstable 
(phase II) and another is meta-stable (phase III). 

• At { = = 0.23125, a GWW type phase transition ^ ES] occurs in 
phase II and the gapped solution becomes a solution with non-uniform 
distribution (Phase IV). 

• At ,^ = ,^2 = 0.237, there is a first order phase transition between the 
phases I and III. Above .^2; the phase III is stable and the phase I is 
meta-stable. 

• At ^ = .^3 = 1/4, phase IV merges into phase I, and the uniform solution 
becomes unstable beyond ^3. 

These are summarised in figure [31 



14 



F Phase II (gapped) 



Phase IV (non-uniform) 

?3 f 




Figure 3: Free energy fl26l) vs C, in the four phases. The gapped and non- 
uniform solutions here are numerically evaluated. Since ^ is a monotonically 
increasing function of temperature (see f l2^ ). the uniform distribution (Phase 
I) is stable at low temperatures and the gapped distribution (Phase III) is 
stable at higher temperature. A first order phase transition between these 
two phases happens at ^2- 



Using (El]), we can read off the critical temperatures corresponding to 
these transition points: 



f3m = mm) 



2Ao 



W 



3(2vre)V3 ^ X 



A 



2\ 2/3- 




(log ] - 1.53 - log ) , m = 0, 1, 2, 3. 



(31) 



Here Lambert's W function is employed again. In the second step we have 
assumed Aq/A ^ 1 and = 0(1). 

As we come down from /3 = oo (go up in temperature), there is a first 
order phase transition at P2 from the centre symmetric phase (TiU = 0) to 
the broken symmetry phase (Trf/ 7^ 0) at an inverse temperature 



1 /A2~ 



(32) 



Several comments are in order here: 



The assumption, used in Section 12. H that e~^^ is small, is correct 
in the parameter region we are interested in. The interesting phase 
structures appear in the regime ^ ~ 0(1). Thus since A/Aq <^ 1 
-A/3 _ g-Ao/3 ^ I f^Qj^ Therefore the terms in 

2, 3, ... are suppressed. 



(see (Uni)), e 
involving U"', n 
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(b) The Euclidean model (E]) is symmetric under the interchange of {t, /3) <H- 
[x, L). Hence, similarly to fl32]) . we can deduce a phase transition in L 
from the TrV^ = phase to TiV ^ (at large enough /3) at a critical 
length Lcr = P er- 
ic) The existence of a finite Lcr above confirms that the transition which 
we found at L — oo and /3 = Per between TtU = and TiU ^ indeed 
happens in the TrV^ = phase. Therefore the expression for Per is valid 
even at finite L as long as Trl^ = 0, since large volume independence 
[HI [im ensures that gauge invariant quantities like the free energy and 
vev of Wilson loop operators do not depend on L in the TiV = phase. 
Thus, the correct definition of 'large L' in this section is 

L ^ Lcr = Per: (33) 

which ensures that we are in the TrV^ = phase. Per is defined in (l32l) . 

(d) By considering the interchange P L, we can claim that, if there is 
no direct transition from TiU = TtV = phase to TrU ^ TrV^ ^ 
phase, the two transition line P = Per and L = Lcr nieet at P = L = Per- 
See figure HlR 

(e) As we discuss in appendix [Cl the non-zero mass in ([3]) does not change 
the qualitative nature of the phase structure. 



3 Phase transitions at small L 

In this section, we discuss the phase structure for small L [L <^ Lcr)- In this 
case we can dimensionally reduce the theory (see footnote H]) to obtain the 
action (l42l) with d = 1. Hence we can use the analysis in [6] where the 
phase structure has been studied by using the 1/D expansion. The phases 
are characterized by the eigenvalue density (127|) . Here we summarise the 
results of [6], 

• P > Pel- The stable solution is given by m„ = (n > 1). The eigenval- 
ues of Aq are distributed uniformly. 

15 Note that a transition line between TiU = TrV = phase and TrU ^ TrV^ ^ 
phase, if it exists (e.g. as in the third circled option in figure |4]), could depend on L and 
/3. However the gravity analysis in Section [5] suggests that there is no such transition in 
our model, consistent with the 'cascade picture' reviewed in Section [Hj In this case, as 
represented by the second joining option in figure 21 can be relaxed to i > Lcr- 

i^In |S] we had considered massless adjoint scalars. We generalize the results to non-zero 
mass in Appendix [Cl 
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• Pel > P > Pc2'- The stable solution is given by ui 7^ 0, u„ = (n > 2). 
The eigenvalue distribution is non-uniform and gapless. 

• Pc2 > P- The stable solution is given by n„ 7^ (n > 1). The eigenvalue 
distribution is gapped. 

• The phase transition at /3 = Pd is of second order and the transition 
at P = Pc2 is a third order (GWW type) transition. 

The critical temperatures are calculated up to 0(1 /D) in [6] as 
PeAr = ^ogD{l + -{^-^] 1, (34) 

Pc2~X\ 

203V5\ , ~ 1127^5 85051 

— \ogD — + 

480 J ^ 1800 76800 

(35) 

where D = D + 1 and Ai = {g'^N^D + 1)/L. In the P-L plane the transition 
lines appear as curves P cx L^/^ passing through the origin. Since our analysis 
is valid only for L <^ Lcr, we should trust these transition lines only in that 
region, as we have depicted in figure HI By using the /3 -H- L reflection 
symmetry, we can also infer phase transition lines for P <^ P^r described by 
L oc as shown in figure |H 

Contrary to the case of large L, an intermediate non-uniform phase ex- 
ists at small L. A similar feature in the context of higher order confine- 
ment /deconfinement type phase transitions has been seen in [20] . 

We should mention that considerations in this section are valid up to 
A' < Xmax where Xmax = L/P^ for P < P^r, and Xmax = P/L^ for L < L^r 
[B], which can be large close to the origin (See footnote [M] also.). We will 
come back to this point in Section |5l 




4 Phases of 2d gauge theory on 

In the last two sections, we have studied confinement /deconfinement type 
transitions in the model ([3]) for large and small values of the spatial size L. 

^''Eqns. (|34 p . (|55|) arc calculated for m — Q. The massive case is discussed in appendix 
O Note that the mass from the KK modes for d = 1 is proportional to VAiL as in (|4T|) . 
Thus the mass correction is small for small L. 
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We have found that the nature of the transition depends on L. We can 
summarise these results in figure HI where we supplemented our calculations 
with the reflection symmetry /3 o L of the model. 




Figure 4: Phase structure of the 2d gauge theory at weak coupling (defined 
below). There are essentially 4 phases characterized by non-zero values of 
various Wilson lines. The inner region, with both Wilson lines non-zero, 
includes 2 additional phases in which the eigenvalue distribution is gapless 
but non- uniform. The orders of the phase transitions (1st, 2nd, 3rd) are 
indicated. Our analysis does not apply to the region enclosed by the dotted 
lines. Possible connections between the phase boundaries across this region 
are suggested in the inset (where boundaries of the intermediate phases are 
omitted for simplicity). A similar diagram is proposed in [3] for the model 
([3]) with large mass for the adjoint scalars (see Section l6T2] for details). As 
we will see in figure the gravity analysis conforms to the second pattern. 
We will see in section I^TT] that the second pattern is also supported by lattice 
studies. 



Weak coupling in the above diagram (figure Hj), for large L (or large /3), is 
defined by A/A^ ^ 1 (see assumption (a) below fITU]) ). In case the 2d gauge 
theory is obtained from a KK reduction, the above notion of weak coupling 
translates to A <^ Mj^j^ (see (|T9l) and footnote [12]) . For small L (or /3), the 
coupling should satisfy A' < f3/L^ (or L/(3^) to validate the additional KK 
reduction to one dimension; as remarked at the end of Section [3l this limit 
on A' can be quite large close to the origin of figure HI 

As mentioned in the Introduction, our model ([3D can be regarded as a 
dimensional reduction of a D + 2 dimensional pure Yang-Mills theory com- 
pactified on a small T^. For the dimensional reduction to work (see footnote 
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Hj), Wj = TtUj {I = d, ■ ■ ■ ,d + D — 1) must be non-zero (which is ensured 
by a sufficiently small size of the D-dimensional torus). Therefore we can 
regard the phase structure in figure [Has a part of the D + 2 dimensional pure 
Yang-Mills theory in the Wi ^ phase. Such a Yang-Mills theory on and 

have been studied in lattice gauge theory and we will compare our results 
with those studies in Section El 

Since our phase structure is derived through the 1/D expansion, it is not 
a priori obvious whether the result should be valid for small D. However, 
at least for small L, the comparison with numerical studies [l6l HZl HH], as 
explained in [6], turns out to be remarkably good even for small D. For exam- 
ple, for D = 2 the 1/D expansion, performed up to an accuracy of 0(1/ 
reproduces numerical results within the expected 25%. Thus we believe that 
the phase structure in the large L region also should be qualitatively correct 
for small D {D>2). 



5 The phase structure from gravity 

In the previous sections, we evaluated the phase structure of the bosonic 
gauge theory ([3]) at weak coupling. In this analysis, it was difficult to figure 
out the phase structure of the middle region, namely where both /3 and L 
have intermediate values. In particular, it was not clear how the various 
phase boundaries in figure S] are connected. 

In this section, we attempt to construct a gravitational dual of our system 
along the lines of Witten's realization of the 4d pure Yang-Mills theory \22\ . 
We consider IIA supergravity on i?^ x (T^ = SjjXSlx SlJ and put N D2 
branes on the T^. The AdS/CFT duality in this context is discussed in [23] 
and more details are provided in Appendix [Dl The dual gauge theory is 3 
dimensional A/" = 8 SU{N) super Yang-Mills on T^, with the identifications 
(t,xi,X2) = {t + f3,xi + L,X2 + L2). In order to complete the definition of 
the theory, we need to choose a boundary condition of the fermions along 
each compact direction. Let us impose an anti-periodic boundary condition 
on the fermions on the X2 cycle. If L2 is sufficiently small (I/L2 ^ A3 = g^N 



^), we can use a dimensional reduction to the two-dimensional torus Sj^x S}^: 
owing to the anti-periodic boundary condition along L2, all fermions would 
acquire a mass proportional to I/L2 and we can ignore theni^. 



i*^This implies I/L2 > 



In fact, even the scalars would acquire a mass at one-loop, as in (35]. However, unlike 
in [35], the scalar mass does not become infinite as L2 — ^ 0. From the 2d gauge theory 
perspective, the scalar mass renormalization due to fermion loops schematically goes as 
ttiy = J '^^P (TO)(p+m) ~ ^^m^ ~ "^^cre we have used a fermion mass m = I/L2 
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One would, thus, expect the gravitational system, for small L2, to describe 
a dual of with D = 8. Unfortunately, however, the gravity solutions 
are not valid in the small L2 region (Lo <g; I/A3) since stringy corrections 



become important (see appendix [DX2l)H In case of L2 ^ V-^s? since the 
fermions are not decoupled, the D2 brane theory will depend on the boundary 
conditions of fermions along the t and Xi directions. There are 4 choices of 
boundary conditions: (AP,AP), (AP,P), (P,AP), (P,P), where P denotes the 
periodic boundary condition and AP denotes the anti-periodic onj^. Phase 
diagrams of the gravity theories for different spin structures are worked out 
in Appendix [Pl and presented in figures [5l [7] and [3 The salient features are: 

(i) The phase structures in the gravity analysis depend on the boundary 
conditions. 

(ii) Only the gravity analysis with (P,P) boundary condition is reliable 
as a prediction for gauge theory through the arguments in appendix ID.4[ 
The phase structure in this case is shown in figure [51 It predicts the second 
joining pattern in figure HI 

(iii) The phase transitions in the gravity description in figureOare Gregory- 
Laflamme (type) transitions [IHl EDI EI] and are expected to be of the first 
order, at least for large L2 [231 E| • 

(iv) The gravity analysis for small /3 and L is not reliable. 

and an uv cut-ofT for the 2d theory A = 1/^2- For a more precise calculation see appendix 
121 Since the scalars remain light compared to the KK scale, we must keep them in the 
Lagrangian as in ([3]) . 

This is a common problem in the construction of holographic duals of non- 
supersymmetric gauge theories. Since the gauge theory coupling constant A3 is greater 
than the KK scale (I/L2) in the region of validity of gravity, the gravity description has 
been likened (cf. [Hj, p. 196-197) to strong coupling lattice gauge theory, the small L2 
limit being regarded as analogous to the continuum limit. Interesting results, including 
the qualitative predictions in [21], and results in AdS/QCD, have been obtained using this 
philosophy. We will use the small L2 extrapolation of our gravity results in this spirit. 
^^The gravity calculation with /3 = L in the (P,P) case has been studied in pSl [5]. 
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Figure 5: Conjectured phase structure of the gauge theory from the gravity 
analysis. We used a large and the (P,P,AP) spin structure of the fermions 
on the three-dimensional torus (with AP boundary condition on the Scherk- 
Schwarz (SS) circle). "YxW is the Polyakov loop operator along the SS circle 
fl77p . The gravity analysis is reliable only in the region above the dotted 
line. The transitions in this diagram are predicted to be first order phase 
transitions. (See appendix ID.4p 



Comparing figures H] and O we can see that both diagrams share some 
common features. Both have four phases in similar parameter regions. In 
particular, the behaviour of the transition lines for large fi and large L is the 
same. The line BC in figure E] is independent of L. This is consistent with 
large volume independence [HI [101 IS2] , since Trl^ = on both sides of BC . 
Similar remarks apply to the line AB as well. 

In the small /3, L region too the two phase diagrams share similarities. 
In figure [5] two phase transition lines emanate from the point D towards 
low values of /3, L. However it is not clear from the gravity analysis how to 
continue towards the origin. On the other hand, in figure S] the region near 
the origin O can be computed reliably and the two (double) lines 0D\ and 
OD2 can be identified as a continuation of the two phase transition lines 
mentioned above. We should note that the phase structure in the small 
/3, L region of figure |4] can be calculated from gauge theory even at strong 
coupling, up to A' < \max where \max = L/f3^ for /3 -C /3cr, and Xmax = 
for L <^ Lcr', as A grows stronger, the calculable region becomes narrower. 

In addition to the above similarities, various details of the phases in figure 
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mandOare also similar. Recall that the phases in the gauge theory are charac- 
terized by three solutions: uniform, non-uniform and gapped. Correspond- 
ingly, three solutions (uniform black string, non-uniform black string and 
localized black hole) play a key role in the discussion of Gregory-Laflamme 
transitions in gravity. For large L, the free energy of these solutions in the 
gauge theory are related as shown in figure [31 A similar relation has been 
found in gravity [53], in the case where the GL transition is of first order. On 
the other hand, if the GL transitions are of higher order, it consists of two 
transitions: a transition between a uniform black string and a non-uniform 
one, and another transition between the non-uniform black string and a lo- 
calized black hole [50]. This is precisely similar to the higher order phase 
transitions in the gauge theory, which we have observed in the small L case. 

An important consequence of the gravity analysis is that we can guess how 
the phase transitions in figure H] are connected. In particular, it indicates that 
there is no direct transition between Trf/ = TrV^ = and Trf/ ^ 0, TiV ^ 
phase. In [5], it was pointed out that this property has also been observed in 
large pure Yang-Mills theories on 4 and 3 dimensional tori in the lattice 
calculations of [21[511[55]. Thus the gravity analysis is in agreement with the 
lattice calculation. More details of the lattice calculation are presented in 
Section [61 

In summary, the full phase diagram of the two dimensional gauge theory 
([3]) may be obtained by combining the results from gauge theory and gravity. 
The result would be given by figure [H with the second joining possibility. [^ 

6 Relation to other works 

In this section, we detail some of the remarks made in the Introduction 
regarding previous works. 

6.1 Comparison with lattice studies 

Large A^ Yang-Mills theories on tori have been studied using lattice methods 
in [21 [511 ESI [1], in (i = 3 and 4 dimensions. Reference [1] contains a nice 
summary of these works. We describe some of the salient features below (see 
also figure [6]) . 

(a) If we start from d = 4 pure Yang Mills theory on a Euclidean torus 

^^The extrapolation involved in this conclusion has additional support from the lattice 
calculations mentioned above. 
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with < L2 < Li < Lq then for all radii large enough the centre 
symmetry Z% (see footnote [3]) is unbroken and all the Wilson loops 
vanish. This phase is called the Oc phase. In this completely unbroken phase, 
the thermodynamics in the large N limit does not depend on any one of the 
lengths illiniE]. 

(b) As L3 is decreased, below a certain value L3 there is a phase transition 
to a new phase where the centre symmetry is broken to and W3, develops 
a non-zero expectation value. The other Wilson loops Wq, Wi and W2 still 
vanish. This phase is called the Ic phase in which there is no dependence 
on the lengths = 0, 1,2 which characterize the directions of unbroken 
centre symmetry. 

(c) If L2 is now decreased, maintaining L2 > L3, a new phase 2^ appears 
below L2 (which is a function of L3) where W2 becomes non-zero. The centre 
symmetry is broken to Z^, with non-zero values of while Wq, Wi still 
vanish. 

(d) Proceeding similarly, a phase 3c is reached when Li is reduced below 
a critical value L1{L2, L-^), and the phase 4^ are reached when, finally, Lq is 
reduced below Lq(Li, L2, L3). 

(e) Thus, d = A pure YM theory (with < L2 < Li < Lq) exhibits a 
cascade of transitions 

Oe ^ Ic ^ 2c ^ 3c ^ 4c. 

(f) In case of d = 3 pure Yang Mills theory (with L2 < Li < Lq) the 
sequence of transitions works similarly, leading to a cascade 

Oc ^ Ic ^ 2c ^ 3c. 

(g) It was found in [M] that the cascade of transitions persists even 
when all radii are the same. E.g. in case of the system in (f) with Li = 
L2 = L3 = L, for high enough L all Wf^ = 0; as L is reduced below a certain 
critical value Lc, only one of the PV^'s picks up a non-zero value [2j, and 
the 3D cubic symmetry group spontaneously breaks down to the symmetry 
group appropriate to a square lattice. 

(h) Generally speaking, it was found in these works that the Wilson lines 
Wfj, change from zero to non-zero one by one; two or more Wilson lines never 
simultaneously change from zero to non-zero values. 

(i) Order of phase transitions^ There is ample evidence that the first 
of the cascade of transitions Zp~^ — Zp~^~^ is first order. In the case 

^■^ Since all directions are equivalent in the Euclidean space, the ordering chosen here is 
arbitrary and all arguments below can be repeated with any other ordering. 

^^We thank Rajamani Narayanan for pointing out some references in this paragraph. 



23 



of Oc — Ic transitions in four- dimensional Yang Mills theories on such 
evidences are presented directly in and indirectly, assuming large vol- 
ume independence, in [571 EH EH]. Evidences for the first order nature of 
the Oc — >■ Ic and Ic — 2c phase transitions have been presented for Yang 
Mills theory on in [60], which indicates that the first two transitions 
Z'^^ Zp^^^^ — > Zp~^~'^ are also first order. Our gauge theory analy- 
sis presents analytic evidence that the first order nature continues till the 
transition Z^j — )■ Z^^ (which is 2c — ?■ 3c in the notation of pure Yang Mills 
theory on T^), whereas the last transition Zn — )■ 1 in which the centre sym- 
metry is completely broken, occurs (in the parameter region of Section [3]E1) 
in two steps through a second and a third order phase transitions. The last 
statement, first derived in [5], is corroborated by the numerical work in [15] . 

Let us compare the above with the phase diagram in figure HI which 
describes phases of the theory ([3]) . Note that the theory ([3]) with D = 2 Ef| is 
precisely the one obtained after the steps (a)-(b) described above, in which 
we reduce a d = A, D = theory (pure YM theory in four dimensions) on 
two small circles of length L2^L^. To make the correspondence more explicit, 
we identify Li = L,Lq = (3, Wi = TrV, Wq = TrU. It is easy to see that 
the phase transitions in the L < P region of figure H] precisely correspond 
to the phase transitions described in (c)-(d) abov^. The top right part of 
figure m (above AB1B2C) represents the 2c phase. The region above the line 
OD1B2A (or its mirror image: the region to the right of the line OD2B2C) 
corresponds to the 3c phase. The enclosed region OD1D2O corresponds to 
the 4c phase. The phase transition across ABi from right to left corresponds 
to the transition 2c — )■ 3c; the phase transition across ODi El from above 
corresponds to 3c 4c (see figure [6]). 

^^In other parameter regions 3c — >■ 4c can be a single first order phase transition [51150]. 
^^We will assume that the essential features of the large D calculation which led to this 
phase diagram remain valid for D = 2. This was certainly true in the d = 1 case discussed 

^''The /3 < L region is the mirror image; it corresponds to the sequence similar to (c)-(d) 
resulting from an ordering Lq < Li. 

^*We are treating the non-uniform phase as part of the 4c phase here. 
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Figure 6: Cascade of phase transitions for pure Yang-Mills theory on a 
(adapted from see the points (a)-(e) above for more details). Our results 
in this paper, for the theory ([3]) with D = 2, describe the indicated phases 
2c, 3c and 4c. The Oc — Ic transition is found to be first order from lattice 
studies in [SB] ; the 2c — ?■ 3c transition is found to be first order from our anal- 
ysis; the 3c — )■ 4c transition for an asymmetric torus is found in our analysis 
to be a double (2nd order -|- 3rd order) transition for an appropriate param- 
eter regime (see figure H]), although, it can be a single first order transition 
at other regimes P, EO] . 



The comments (g)-(h) above have a direct bearing on the possible joining 
pattern in figure HI Out of the three possible joining patters shown in the 
insets, the first and the third patterns allow a direct transition 2c — )■ 4c and 
are, hence, inconsistent with the comment (h). Thus, consistency with the 
lattice results described in this subsection uniquely pick up the second joining 
pattern. As we showed in Section [5|, the same joining pattern is also picked 
up uniquely in figure O through the analysis of the gravity dual. 

A more quantitative comparison of our work with the above lattice studies 
is left for the future. 

We should mention another important lattice study [61] which deals with 
super Yang Mills theories in c? = 2 and is closely related to the work presented 
here and in [B]. In parameter ranges where the two theories coincide, our 
phase diagrams agree (see section 5 of [SI])- See also related numerical works 
about the center symmetry breaking in super Yang- Mills [621 [63] . 

6.2 Comparison with earlier analytical studies with 
massive adjoint scalars 

Reference [3] considers the theory ([3]) in the limit m ^ A2^^^. Our figure H] 
is similar to figure 13 of [3], except that our figure is obtained for any mass 
(including m = 0) where their figure is for the large mass limit. The reason 
for the agreement is the appearance of a dynamical mass A for the adjoint 
scalars in our model, as we have explained above. 
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In figure 14 of [3] an interpolation between small and large radii is pro- 
posed on the basis of some analytical estimates for the intermediate radii. 
Our figure O although similar to this figure, differs in one crucial respect. 
The phase transition line BC in our figure is horizontal throughout, as it 
must be according to large volume independence arguments [9l|T0]. Since 
the line BC is entirely in the Ti V = phase, the transition temperature Pc 
cannot depend on L; hence BC must be horizontal. This property is violated 
by the corresponding line (the intermediate radius segment) of figure 14 in 
[3], which should have been horizontal according to the above argument. 

6.3 Comparison with Yang-Mills theories on different 
topologies 

We have found that the nature (in particular, order) of the confinement / deconfinement 
type transition at Pc in the two dimensional gauge theory 02]) at a fixed radius 
L depends on the value of L. Since this theory can be obtained from a pure 
Yang-Mills theory on a small x S}^ x S^, it is interesting to compare this 
result with Yang-Mills theories on compact spaces with other topologies. We 
present such a comparison in Table I6.dp^l . 



M 


type of phase transition 


small T^x small 


2nd+3rd 


small T^x large 


1st 


small 5*^ 


2nd+3rd 


small 5"^ 


1st 



Table 1: Confinement /deconfinement type transitions in pure Yang-Mills 
theories on S*^ x A^. Here "small S^" and "small T^" refer to sizes small 
enough to ensure (a) that the Zjsf symmetries in the and directions, 
respectively, are broken, and (b) that all the KK modes can be integrated 
out. "Large S^" ensures that the symmetry along the is not broken. 

Because of the difficulty of the analysis of Yang-Mills theory, only weakly 
coupled Yang-Mills theories on and 5*^ [511 ES] have been studied. In these 
cases, all the spatial components of the gauge field have a mass proportional 
to 1/R, where R is the radius of the sphere. These massive gauge fields can 
be integrated out perturbatively if the radius is sufficiently small [RAqcd <^ 
1). Reference [BH ES] derived effective potentials for Aq up to three loop 

^^In order to apply 1/D expansion, we need _D > 1 in the small T^x small case and 
D > 2 in the small x large case. 
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order and found the transition in the case to be first order [M]. On the 
other hand, the transition in the 5*^ case consists of second and third order 
transitions as we found in the small case [65] . Note that the higher order 
transitions for small S"^ is expected to change to a first order transition in a 
strong coupling regime according to lattice studies [66] , 

Thus the nature of the transition depends not only on the size but also 
on the topology of the compact space. It would be interesting to understand 
the origin of these differences. 

7 Conclusions 

In this work, we have computed the phase diagram of two dimensional Yang- 
Mills theory with adjoint scalars ([3]), which can be obtained from a KK re- 
duction of a higher dimensional pure Yang- Mills theory. We treated the case 
of massless adjoint scalars in detail, outlining the generalization to arbitrary 
non-zero mass in Appendix O, and found the phase diagram in figure |H At 
large spatial radius, there is a first order confinement/deconfinement phase 
transition, whereas at small spatial radius, there are two closely spaced phase 
transitions: (a) a second order phase transition from the 'confined' phase to 
a 'non- uniform' phase (non-uniform eigenvalues of the Polyakov line), fol- 
lowed by (b) a third order phase transition from the 'non-uniform' phase to 
a 'gapped' phase. Our calculations, based on the large D method [B], provide 
an analytical derivation of the dependence of the thermodynamic behaviour 
on the size of the spatial box, which is anticipated on the basis of lattice 
studies and gauge/gravity duality. 

We have also considered the phase transitions in the gauge theory from 
the viewpoint of a gravity dual, based on a scaling limit of Scherk-Schwarz 
compactification of a D2 brane on a 3-torus. Although there is no overlapping 
region of validity of the gauge theory and gravity descriptions, the analysis 
of the gravity dual leads us to conjecture a certain specific completion of the 
phase diagram in the gauge theory, as in figure [51 In performing this analysis, 
we encountered an inherent problem with the holographic analysis, namely 
a dependence of the physics on the fermion boundary conditions, which was 
absent in the gauge theory description (see Appendix [D|) . Indeed this prob- 
lem is related to a more general problem in the holographic description of 
QCD [67]. We discuss this problem further in |68j . 

We matched our findings from gauge theory regarding the order of various 
phase transitions with those from a gravity analysis in Section \5\ and with 
those from lattice studies in Section [61 

Note that the method of integrating out the adjoint scalars using a. 1/D 
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expansion works equally well in higher dimensional {d > 3) gauge theories 
fl42p . leading to an effective action for the gauge field as shown in (jUD- 
However it is difficult to evaluate the dynamics of this model because of the 
existence of dynamical gluons. This is a crucial difference from the lower 
dimensional cases {d = 0,1,2). In addition, the d dimensional model (142|) 
typically appears through a KK reduction of a d + D dimensional (super) 
Yang Mills theory, but for d > 3 the mass of the adjoint scalars induced from 
loops of KK modes is large (see appendix hence the contribution of the 
adjoint scalars may be not relevant for d > 3. 

The investigations in the present paper were partly motivated by a desire 
to understand a gauge theory dual to a dynamical Gregory-Lafiamme transi- 
tion. The considerations in this paper provide a step towards understanding 
this issue; details of this will appear elsewhere [27] . 
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A Mass for from one-loop of the KK modes 



If we consider a d + D dimensional Yang-Mills on T'^+^ and consider a 
dimensional reduction by taking the radii of to be small, we will classically 
obtain a d dimensional gauge theory coupled to D massless adjoint scalars. 
However, if we consider quantum effects, the action would be modified. One 
of the relevant corrections is that the adjoint scalars would acquire mass as 
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in In this appendix, we evaluate the mass at a one-loop level PI. 

Starting from the d + D dimensional action ([T]), we can derive a one- loop 
effective action for the constant diagonal components of A'^ = (a^*, ■ ■ ■ , a^) 
by integrating out all the other modes [H [691 E] , 

/d+D~l 

seff=- \ n ^ 



' d + D-2T{{d + D)/2) 



/x=0 



7l{d+D)/2 



X 



exp {i Y.U kuLyi^a'i - aj)) 



i,j {fc^}^{o} IZ^i. '^i.-^i/J 



Here the sum J2{k^}^{o} includes all integers except kQ = ■ ■ ■ = kd+n-i = 0. 
Let us now take (/x = 0, ■ ■ ■ ,d — 1) large and Lj {I = d,- ■ ■ ,d + D — 1) 
small and derive a low energy effective theory by using this expression. Gauge 
invariance implies that the effective action in this situation will be given by 

= Iy^ f/'^' (i^- + 1 1 (D.Y^T - E '-fr^ y'w'^ y 

_ f^'ri f'^, . \ d + D-2Tiid + D)/2) ^ \Tre^^^^.^^^-^y' 

JO i=i-J^ \i=i J {fe/}7^{0} \I^J'^J''J) 

(37) 

Here we have rewritten A^+i^i = g^Yj and L^+i^i = h- Qd and g'^ are the 
same as g and of ([3]); they satisfy fi'^+£)/ Il/Li h = Qd = 9d ^ physical 
scale /i > 1/// (c/. (USD). 

If // are small and the long string modes are suppressed (see footnote H]), 
we can treat perturbatively. Then we can expand the exponentials in f l57|) 
and obtain a quadratic term in Y as 



\/=i / {ki}^{ 



km TriY'f 



I'V U2]2\{d+D)/2 2 

}^{o} \2^J '^J''J) 

(38) 

Other interaction terms from the exponential would be suppressed by small 
//. If we take all the li to have a common value Ikk = ^/M^k, we obtain a 
mass for Y^ as 

,'NM'-Hd\D o /((^ + ^)/2) ^ I 1 

(39) 
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We do not use the large D limit in this appendix. 



29 



This sum would diverge for d < 2 and hence needs to be regulated e.g. by 
using a prescription (/c)^ = kf + ■ ■ ■ + kjj < {Kuv / Mkr)"^ , where Kjjv is a 
cut-off scale Ell- This would imply a non-trivial RG flow of the mass. Let us 
first consider the case of c? = 2. For Kuv ^ Mkk, the regulated sum (139|) 
approximately gives m'^{Auv) = A g^N \og{Kuv/MKK) where A is a numerical 
constant. Let us choose a renormalization scale > 0{Mkk)- We can then 
define a renormalized mass at the scale fi as 

= m2(A„,) - A gjN log(A„,//i) = A'Aa, (40) 

where A' = A\og{fi/MxK)i ^2 = 92^- Note that the running of the mass 
will stop below /i < Mkk- 

For d = 1, a similar analysis gives 

^ X^/Mkk. (41) 

For d > 3, the sum in (jSHD is convergent, leading to m ~ X^J"^ M^^^^^"^ 
which is much larger than the typical QCD scale, e.g. for d = 3 the QCD 
scale is ©(As) whereas the adjoint mass is 0{\]/'^mII'^). Thus if M^k is 
large, the adjoint scalars will not contribute to the QCD dynamics. It means 
that only the d dimensional gauge field dominates the dynamic^. 

So far we have considered the mass correction from the KK modes in 
pure Yang-Mills theory. We can extend this calculation to the KK reduction 
of supersymmetric theories with a Scherk-Schwarz compactification e.g. to 
the derivation of ([3]) with D = 8 from the D2 theory. In this case, we need 
to evaluate the contribution of loops of adjoint scalars as well as fermions. 
However, it can be shown that we obtain a similar mass ~ 0{Xl/'^ M^^'^^^'^) 
even in this case [5]. 



B Derivation of effective potential (1111) 



In this appendix, we show the derivation of the effective action flTT]) from 
ffTU]) with the assumptions (a) and (b) below fITU]) . We can evaluate (fTU]) in 



An ultraviolet cutoff typically breaks gauge symmetry. The calculation discussed here 
can be repeated avoiding such problems, by using dimensional regularization [7]. 

■^^For d = 2, the mass of the adjoint scalar from the KK modes is finite (0(A2)) but 
the dynamical mass A is (logarithmically) larger than A2 (see HH)). Hence one may 
naively think that the adjoint scalar would be irrelevant as for d > 3. However this is 
not correct [3] , since the phase structure of the 2d pure Yang- Mills on is trivial and is 
always confined. Therefore the contribution of the (logarithmically) heavy adjoint scalars 
is important in the 2 dimensional gauge theory 
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a general d dimensional gauge theory, described by the action 

/,J / 

This model can be identified with f p7|) if we take Lj small and choose the mass 
and couplings appropriately. We will first discuss the general rf- dimensional 
case and apply the results to c? = 2 later. 

Let us first set m = 0. Then, through a similar calculation as in Section 
12. H we obtain a generalization of ( ITU]) (the calculation closely follows [3]). In 
this section, we will use the more general notation {Lq,Li) for {f3,L). 



6S[A, A] =^ log det {-Dl + A^) = ^Tr log f E + A, 



{n^} 

D Lq - ■ ■ Ld-i 



2 

+ A^ 



2 (2ttY 



Tra., {L,})e^J^.^-^''^^ (43) 



where P{fc^} is 
P|,,4(A,{L4) 



{n,} Va.=0 V / / 



/no+27r/Lo pnd-i+2Tr / L^-i 

dao--- dad-i log E + A' ) 



da,--- daa-i log E «S + A2 e-* E, k,L,a, 

OO J —OO \ , n / 



(44) 



Let us now evaluate P{k^} in the {k^} = {0} and {k^} ^ {0} cases separately. 
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In the {k^} 7^ {0} case, P{fc^}(A, {L^}) becomes, 
P|,^}(A,{L4) 



dar 



-00 
00 



'd~l 



dao ■ 



dad-i log [ ^ aj + A^ j c-^^m'^m^'m'^. 

At=0 



-00 ^ —00 

Jo 



dad-i lim 

e-5.0 



-loge - 7 



ad/2+1 



2 (27r A) 



d/2 



d/2 



^1 



2(27rA) 



d/2 



d+l 



TT 

2A 



exp 



(45) 



where is the modified Bessel function of the second kind and we have 
used Ka{z) = \/ Ti I2ze~^ + ■ ■ ■ for large z in the last equation. 
For {k^} = {0}, we find 

/oo POO ( 
dao--- dad-i\ogiJ24 + ^ 
'OO J —00 \ n 



/x=0 



daQda"^-^ log (a^ + A^) 



(46) 



where ^Id = 27i'^^'^ /T{d/2) is the surface area of the d dimensional unit sphere 
and A is a cut off. 

Finally we turn on the mass term. In this case, we can obtain the results 
by replacing A — >■ a/A^ + in (H5|) and (146|) . Note that the assumption (a) 
and (b) below ( JTOj) should be modified accordingly. 

B.l Effective action for d = 2 

We now consider the special case of d = 2. In this case, we can evaluate P{o} 
as 

.A 

P{o} (A, 13, L) =2tx da a log [a^ + A^) 
Jo 

= - ttA^ + TT (A^ + A') log (A^ + A^) - ttA^ log A^. (47) 
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By using this result and fHSl) to fH51) . we obtain the effective action flTT]) . Note 
that Tradje*^'^'^^^"'"'-^^^^ in becomes |Tr (f/^'V') |^ by using the assumption 
(b). 



B.2 Effective potential for higher dimensional models 

It is easy to generahze the derivation of the effective potential f l2U]) for the two 
dimensional gauge theory in Section 12.11 to the d dimensional gauge theory 
2]) for large o By using the results in the previous section, we obtain 



11=0 



1 / , d+i ^ 



Po {\l^l + m\ {L,} 



Here Aq is defined as a solution of the saddle point equation 



A2 
2A' 



1^. 



da- 



-Tre*'^''^'' 



(4^ 



(49) 



{2'kY Jo a2 + A2 + m2' 
Note that we can always find a unique positive solution Aq from this equation 



C The phase structure of massive model. 

In this appendix, we study the two dimensional gauge theory with a 
mass term for the adjoint scalars. As we mentioned in the introduction and 
elaborated in appendix |Al such a mass term generically arises from KK loops. 
We discuss here how the results of the massless case in Section [2] and Section 
E] are modified. We will show that the mass does not change the qualitative 
nature of the phase structure. 

C.l Large L case 

By using the results in IB.2[ we generalize the effective action for TtU (12 3 p to 
the massive case. The resulting effective action is again given by (12 3 p with 

■^■^Tlic 1/D expansion in the d dimensional gauge theory (|^^ in a high temperature 
region is considered in [37] also. 
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different values of ^ and C: 



27rA2/33 



C(A',Ao) 



57r 



1 



A' 



57r 



27rA2 



(50) 
(51) 



Tlie dynamical mass Aq is a solution of 

A' = - 



27rA2 



A2 



(52) 



For large A, Aq becomes 



Ar 



/ A' ^ |^27rA2 



27r V A' 



-e A' 



(53) 



Therefore we can use the same analysis as in Section [2^21 and obtain the same 
phase structure with the following modified transition temperatures 



(3„ 



W 



3(2vre)V3 



A^ 

^0 



m 



A 



2/3 



(54) 



Although the mass changes the explicit values of the transition tempera- 
tures and some other physical quantities, the qualitative nature of the phase 
structure is not modified, as we have mentioned. 

Note that fl52l) implies that Aq becomes smaller as m increases for fixed 
A' and A. Therefore, for heavy mass, we can ignore the dynamical mass Aq 
compared to m and our calculations reproduce the heavy mass QCD results 
in [HE]. 



C.2 Small L case 

If L is small enough in dH]) , we can integrate out all the non-zero momentum 
modes in the L direction (see footnote Hj) and obtain a one dimensional model 
(H2|) with d = 1. In this case, the mass of the adjoint scalars induced at one 
loop from the KK modes is proportional to (AiL)^/^ (see (14T|) ). Hence we 
can ignore it for small enough L and the results in [H] shown in Section [3] 
would still be valid. Although the contribution from the mass would be 
small, it may be valuable to confirm that the results in [6] are not modified 
qualitatively. 
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Starting from fH21) with d = 1, we obtain an effective potential tlirougli a 
similar calculation as in 



oAi 1 V-'^ / 

^ n=l ^ ' 



n 
(55) 



Here the third term is a contribution of the Vandermonde determinant [B]. 
Then the saddle point equation for becomes 

+ (-=L=e-"^v^)Kr (56) 



The solution of this equation at low temperatures (e-/3V^^+™^ ^ 1) is 



where 



2 

m 



g (/(m) + /(m)-i-l) 



/(m) =^T7^ - 2m' + y/27A?(27A?-4m6)j ^ . (58) 

Note that, although /(m) is a complex for > 27A2/4, Aq is always real 
and positive. By substituting this solution to (155|) . we obtain an effective 
action for 

Seff{{un})/DN' = + + + a\u,\' + b\m\' + ..., (59) 

8A1 ^ 

where 

1 _ b = ^ -2/3v/XlW (QQ) 

D J' 3Al + 2m^ ■ ^ ' 

One important fact is that b is always positive. It has been shown that, in 
this case, the confinement/deconfinement type transition always consists of 



■^"^Contrary to the d> 2 case, the condition (a) and (b) below ((TU)) are not required for a 
derivation of the effective action in the d — 1 case. Thus the 1/D analysis would be valid 
as long as the effective dimensionless coupling does not scale with D. 
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the two transitions (2nd+3rd) as in the massless case P, [SO]- These crit- 
ical temperatures are given by the solutions of a = and a /{2b) = —1/4 
respectively, and evaluated as 



Pel 



1 





+ m- 



/3c2 =/3cl 



Ai logD 
2D 3 A2 + 2m2 ' 



(61) 



Therefore as in the large L case, the qualitative nature of the phase transition 
is not modified by the mass ternjfl. 

Although we have evaluated only the leading order of the 1/D expansion 
in this section, the result does not change even if we include the next order. 

D Phase structure of D2 branes on a 3-torus 

In this appendix, we consider the gravitational system of Section [S] in detail. 
We first review some generalities for Dp branes with arbitrary p. 

Note that our results based on the 1/D expansion disagree with the results in [5]. 
Reference [3] studied the same model (|55l) by using a large mass approximation, which 
is supposed to be valid if Xi/m^ <C 1, up to three- loop order and concluded that the 
confinement /deconfinement transition in this model would be a single first order transition. 
The difference would presumably arise from the fact that the 1/D expansion employed 
here evaluates the model in a non-trivial vacuum characterized by the non-zero A, whereas 
the large mass analysis of [3| is performed as a perturbation around the trivial vacuum. 
Numerical studies analogous to [1^ but performed for massive adjoint scalars [JT] should 
be able to provide further insight into this issue. 
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D.l Dj9 branes wrapped on p+ 1-Torus 
D.1.1 The solutions 

The geometry of a black Dp brane on a p-torus Xj e (0, Li), (i = 1, ■ ■ ■ ,p) in 
the Maldacena limit (assuming Euclidean time t e (0, /3) ) [23] is given by 



a 



F{u) ( f{u)dt^ + dxidxi) j + 



dv? 



i=l 



F{u)f{u) 



G{u)dn 



2 

8-p 



(27r)^-PA 

AT 



\/ dp\p+i 

27-2p^{9-3p)/2p 



G{u) = y%\:;[u^P-^^/\ f{u) 



Uo\7-P 

u 



7 — p 



A 



p+i 



(62) 



where gp_^_i is the p + 1 dimensional YM coupling, which, in the bulk theory, 
can be regarded as specifying a boundary condition for the dilaton field. The 
dimensionless YM coupling, defined by 



9!ff = {9liN)u^-' = (e^N) 



2/(7-p) 



(63) 



is given directly in terms of the dilaton; its dependence on u for p ^ 3 reflects 
the running of the gauge coupling. The scalar curvature is given by 



a'R = 1/g. 



eff- 



(64) 



Since time is Euclidean, u G (uo,oo). The smoothness condition at n = uq 
relates the inverse temperature /3 to uq as follows: 



2n~ 7-p 



(65) 



The classical action of the black Dp brane is [211 [231 E] 



p-3 



S/N' =CpX;-'{Li ...Lp/3i-r^ + H,,^{U) 



_2{7-p2 



Cn 



5 — p 



2ii-2p7r(i3-3p)/2r((9 _ /2)af'-P'>/^'-P^ ' 
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2a„ 



dp = 2^^2P^(9- 



-3p)/2-p ( J^^ ^ 



2(7-p)/(5-p) 



u 



7-p 



(66) 



P+i 
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which is evaluated by expressing the on-shell action as a regularized integral 
of ^/ge~'^'^R over the range uq < u < U. 

In the following we will also be interested in AdS solitons [22l [70] which 
can be obtained from (162|) by moving the coefficient f{u) from dt"^ to one of 



the Xj's , say to dx"^ 



ds^ = a' 



p-i 



i=l 



F{u) I f{u)dxl + dt^ + y^ dxidxi) 1 + , + G{u)dn^ 

> ^ ^-^ I r{u)j{u) 



8-p 

(67) 



The smoothness condition at n = uq now gives a condition analogous to f l65|) 
where /3 is replaced by Lp. Thus this solution has a contractible Xp-cycle 
(which wraps around a so-called 'cigar' geometry on the {xp, u) plane) along 
which the fermions must obey the anti-periodic (AP) boundary condition. 

The regularized classical action evaluated on such a classical configuration 
is given by 

S/N'' =Cp\^,L, ■■■Lp/3 (^-L;^ + i/rcg(f/)) , (68) 

where the notations are the same as before. 

In a toroidal, Euclidean, spacetime, the time direction is on a similar 
footing as any other direction. Thus, the difference between the black brane 
and the solitonic solution is only in the labelling of the contractible cycle 
(location of the 'cigar' geometry). Hence we will sometimes refer to both 
black branes as well as AdS solitons as just Dp solutions. In the following 
sections, we will consider different Dp solutions (with p = 0, 1, 2) which wrap 
on (are localized along) various cycles and, in order to distinguish them, we 
will use the following notation: 

DpLoiLi,...,Lp) denotes a Dp solution which (i) has a contractible Lq cycle 
(that winds around the 'cigar') and (ii) wraps on the Li,...,Lp cycles; this 
is a black brane. Similarly Dpl^(^Lo,Li---,Lp-i) denotes an AdS soliton in which 
the roles of t and x^ are flipped, as in fl67|) . 



D.1.2 Validity of supergravity 

The solutions described in the previous section are leading order supergravity 
solutions. When we consider a black Dp brane solution (l62i) . the gravity 
analysis is reliable, if the parameters satisfy the following conditions [51] : 

1. The typical length scale of the black Dp brane near the horizon (see, 

1 /2 

e.g. (IMl) ) is given by I = (a'y^ dp Ap+in[f ^'^^^^ j . In order to sup- 
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press stringy excitations, we should satisfy / ^ \fa' . From this 
condition is equivalent to 

Ap+i/3=^"^>l (p<5). (69) 

2. The mass of the winding mode along an Li cycle is given by M^i = 

(^a'uQ ^''^^/ ^/dp\p+i^ Li/a'. In order to suppress the winding mode, 
we must have M^jil ^ 1. This condition gives 

AJ/.'a'^»/3 ip<5). (70) 

If this condition is violated and if the fermion on the brane satis- 
fies the periodic (P) boundary condition along the Lj cycle, we can 
perform a T-duality along this direction and reassess the validity of 
supergravity in the dual framj^. After the T-duality, the black Dp 
solution becomes a smeared black D{p — 1) brane solution, which is 
composed of uniformly distributed D(p — 1) branes on the dual Li cycle 
[231 EH Eg. Then the condition ([70]) is replaced by M^il > 1, where 

— 1 /2 

M^i = (^a'uQ^^^^"^ / dpAp+i j (27r)^/Lj is the mass of the winding 
mode on the dual cycle. This condition gives 

U < (3. (71) 

Note that the first condition (169 p does not change under the T-duality 
(the value of the classical action is also invariant). 

If, instead of a black Dp brane, we consider a solitonic solution which is 
obtained from the black brane by flipping t -H- Xi, then the conditions for 
the validity of supergravity are simply obtained by replacing /3 and Li in the 
above conditions. 

We will discuss these criteria below in some detail in the parameter regime 
of interest. 



D.2 Phase transitions of the Dp solutions 

Using the on-shell classical actions ( 166|) and ( 168|) . we can determine various 
phase transitions [25| [26| |5] , as we will show now. 

■^^If the fermion satisfies an anti-periodic (AP) boundary condition along the Li cycle, 
the theory is mapped to a type theory through the T-duality. Then the bulk theory 
involves a tachyon and how the holographic description of gauge theory works in this frame 
is unclear. We thank Shiraz Minwalla for pointing this out. 
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D.2.1 GL transition 



The Gregory-Laflamme (GL) transition [l9l [SQl [5T] was originally found in 
the context of black rings in D > 5 which were found to be unstable unless 
wrapped on a sufficiently small circle. In the context of black p-branes, the 
GL instability shows up as follows [231 ED E2j- Suppose a black p-brane 
is wrapped on a circle of length Li, along which the fermion satisfies the 
periodic boundary condition. If Li is so small that it violates ([70]), we need 
the T-dualized description in terms of a uniformly smeared black {p — 1)- 
brane on the dual circle of length L[ = (27r)^/Li. If L[ is large enough, the 
smeared black (p— l)-brane undergoes a GL transition leading to a black p—1 
brane localized on the dual cycle. The transition can be studied dynamically, 
as well as thermodynamically. To study the latter, let us interpret f l^B]) 
without the regulator term (see Eq. (9) of [21]) as the Euclidean action 
(above extremality) of the smeared black {p — l)-brane : 

Sp/N^ = -CpXfiL, ■ ■ ■ (72) 

Here the action is evaluated in the Dp-frame (recall that the action is invari- 
ant under T-duality). 

The value of the Euclidean action (above extremality) for the localized 
black p—1 brane solution is approximately (see below) given b\l^ 

Vi/^^' = -Cp„iA|^L2 ■ ■ ■ (73) 

For small enough Li (large enough L[), this is smaller than (172]) . The tran- 
sition between the uniform and localized p—1 branes happens when fl72p 
equals fl73]) : 



(3 fa 



(5-p)(6-p)/4 



Ap+iLr^ (74) 



Here we have used Xp = Ap+i/Li. In (175]) one has used the approximation 
that the horizon size is much smaller than 1/Lp. Since this is strictly not 
true near the phase transition point (where the horizon size is of the order 
of 1/Li), the estimate fITi]) is approximate. See [5] for some details of this 
approximation. 

Note that there are several arguments that this transition would be first 
order [231 EH [5]. 



^^Improvements to this approximation are discussed in |71j 
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D.2.2 GL type transition in the soliton sector 



For bosonic theories in Euclidean spacetimes, the (3 Lp interchange is a 
symmetry, provided all other radii are left unchanged. This is a symmetry of 
fermionic theories as well, provided the spin structures along t and respect 
this symmetry. Thus, if there is a GL transition given by fl74p . there must 
be a GL type transition between a uniformly distributed solitonic (p — 1) 
brane to a localized solitonic {p — 1) brane when 1/ Li becomes too large. 
The transition is given by fl74|) with a /3 -H- Lp interchange: 



D.2.3 The Scherk-Schwarz (SS) transition 

This is a transition between the black p brane configuration flB^ and its 
solitonic counterpart ( 167|) . If we use the same large U regulator for both the 
solutions, then ( 166|) and ( 168|) can easily be compared. The regulator term is 
the same in both the actions and can be ignored while comparing the two. 
Equating the two Euclidean actions, we get the transition temperature 



Similarly, we can consider another solitonic solution by replacing Lp L^, 
if Lfc is also an AP circle. Then further transitions will happen at /3 = 
and Lp = Lk. 

Note that this transition is also expected to be first order [26] . 

D.3 D2 branes for various spin structures: generalities 

In this section, we apply the general properties of the Dp solutions stud- 
ied in the previous sections to the D2 brane on the 3-torus: (t, xi,X2) = 
+ xi + Li,X2 + L2) El and make a prediction for the phase structure of the 
2 dimensional gauge theory ([3]). As mentioned in Section [5l we fix AP bound- 
ary condition for fermions along the X2 circle. Then, there are 4 choices of 
boundary conditions for the fermions along the t and Xi directions: (AP,AP), 
(AP,P), (P,AP), (P,P), where P denotes the periodic boundary condition and 
AP denotes the anti-periodic one. We will evaluate the phase structure of 
the gravitational system with these 4 boundary conditions. (Since the result 
in the (P,AP) case can be obtained from (AP,P) by exchanging /3 ^ Li, we 
will show the results in the (AP,P) case only.) 

■^^The notation Li in this section represents what is called L in the main text. 




(75) 



l3 = Lp. 



(76) 
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The order parameters of this theory are the Wilson loop operators winding 
around each cycle: 



If the gravity solution has a contractible cycle (i.e. it wraps around a 'cigar'), 
the expectation value of the corresponding Wilson loop operator is non-zero. 
If the solution is localized on a cycle, then also the expectation value of the 
corresponding Wilson loop operator is non-zero. However, if the solution 
wraps around a non-contractible cycle, the expectation value of the corre- 
sponding Wilson loop operator vanishes [23] . 

In order to derive the phase structure corresponding to figure HI we will 
evaluate the phase structure of supergravity for each boundary condition by 
changing /3 and Li with a fixed L2, which is related to a cut off scale of the 
2 dimensional gauge theory ([3]) (see footnote 120]) . 

From now on, we use units such that A3 = 1. 

D.3.1 D2 on (AP,AP,AP) torus 

We consider the phase structure of D2 branes on a 3-torus with (AP,AP,AP) 
boundary conditions. In this case, three solutions appear: D2p(^L^^L^), D2i^(^p^i^) 
and D2L^(^i3^L-^y (Recall the notation at the end of section lD.l.ll ) In this case, 
the theory does not have the P circle and only the SS transition (1761) happens. 
The phase structure is shown in figure [71 

As we have argued in appendix ID. 1.21 we need to check the validity of 
the gravity solutions. Let us consider the solitonic solution D2i^(^p^L^y From 
(169|) . L2 ^ 1 is required (in units where A3 = 1) to suppress the stringy 
excitations around the tip of the cigar. We need to check the condition 
related to the winding modes also. In order to suppress the winding modes, 
L\^'^ ^ 13'^/'^ ^ L2 are required from fl70l) . The phase boundary AB and BC 
are given by Li = L2 and f3 = L2, and these conditions are satisfied on these 
boundaries in the large L2 case. Thus the D2i2(/3,-Li) phase is always reliable, 
if L2 is large. 

Next we consider the black brane solution D2p(^i^^i^y From (l69l) . /3 ^ 1 
is required. Thus this solution is not reliable when /3 ~ 0(1). We can 
see that the condition (1701) for the winding modes is satisfied in the region 
surrounded by the phase boundaries and (3^1. Therefore the D2p(^L^ L^-) 
solution is reliable in the region indicated in figure [71 Similarly the solitonic 
solution L'2j,j(^ 2.2) is reliable in the region indicated in figure [71 




(77) 
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Tr[7 = 
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TrVK = 



A 



Trf/ = 
TrV = 
Triy ^ 




B 



'c 



7 



Trf/ 7^ 
TrV = 
Trl-y = 



o 



Figure 7: Phase structure of the D2 brane on with (AP,AP,AP) boundary 
condition for large L2. The gravity analysis is valid above the dotted line. 



Summing up these tests for the validity of the gravity analysis, the phase 
structure is reliable if L2 is large and /9 and Li are above the dotted line in 
figure [71 This is, of course, a problem, since we are interested in the results in 
the L2 — limit as we mentioned in Section O We will discuss this problem 
in appendix ID.4[ 

D.3.2 D2 on (AP,P,AP) torus 

In this case, 4 solutions appear: Dl^ij^^^j^^-^.DIi^^i^^^i^^-^.Dli^^i^^-^ and D1^(L2). 
The phase structure for large is shown in figure [HI The phase boundaries 
are given as 



Note that AB is the GL type transition ( !75|l . BO is the GL transition ([74]) 
and EC is the SS transition (1751) . 

Through similar tests for the validity of gravity as before, we find that 
the gravity analysis is valid only in the region above the dotted line in figure 
m (Again the gravity analysis in a small L2 is invalid.) 




(78) 
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Figure 8: Phase structure of the D2 brane on with (AP,P,AP) boundary 
condition for large L2. The gravity analysis is reliable only in the region 
above the dotted line. 



D.3.3 D2 on (P,P,AP) torus 

In this case, 4 solutions appear: i522,2(/3,Li), -D1l2(^), -D1/,2(Li) and DO^^. The 
phase structure for a large L2 is shown in figure O The phase boundaries are 
given by 

AB:L,= (I) ' Lf , BD:fi = L„ : /3 = (|) Lf . 

(79) 

These are GL type transitions fl75|] . Other lines can be obtained by /3 o Li. 
The gravity analysis is valid only in the region above the dotted line in figure 
[5] for large L2. 

We will adopt the phase structure in this boundary conclusion as a pre- 
diction for the gauge theory. The reason will be explained in the next section. 

D.4 Conjectured phase diagram of 2D bosonic gauge 
theory 

In the preceding subsections, we have obtained various phase structures from 
gravity for different spin structures {i.e., for different fermion boundary con- 
ditions). Since these results are valid only for large L2, we need to extrapolate 
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them to small L2, where the bosonic gauge theory should appear (see foot- 
notes [2D] and [22]) . 

From the viewpoint of the two-dimensional bosonic gauge theory, a holo- 
graphic correspondence with the various gravity phase diagrams described in 
this section is problematic since the latter have a dependence on fermionic 
boundary conditions, while no such dependence obviously exists for the two- 
dimensional gauge theory (such dependences, however, exist for the three- 
dimensional SYM theory, which is more directly related to the gravity de- 
scription). Indeed one such problem was pointed out in [67]. This argument 
is further developed in [68]. The arguments presented in these papers indi- 
cate that the phase transition in the bosonic gauge theory cannot be regarded 
as a continuation of the SS transition of gravity. Thus, in order to have a 
smooth continuation of phase boundaries between the holographic descrip- 
tion and the two-dimensional gauge theory, we should avoid choosing the 
(AP,AP,AP) and (AP,P,AP) boundary conditions, in which the SS transi- 
tion appear. For this reason, we choose the (P,P,AP) case to read off the 
predicted phase structure from gravity. 

Note that all the transitions in the (P,P,AP) case are of the GL-type, and 
are supposed to be first order phase transitions [23], [51], [5] . In this case, in 
addition to the uniformly distributed solitonic {p — 1) brane and localized 
solitonic (p — 1) brane discussed in ID. 2. 2] non-uniformly distributed solitonic 
(p — 1) brane, which is always unstable, appears. These three solutions would 
correspond to the uniform, non-uniform and gapped distribution in the gauge 
theory (see figure ([2])). Indeed, the free energies of these gravity solutions 
are expected to satisfy a similar relation to those in the gauge theory shown 
in figure [3] through numerical study in general relativity [53] . This fact also 
supports our prediction from the gravity analysis. 
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